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1. Introduction 

In a previous paper |J the flow of N = 2 SUSY pure SU{2) Yang-Mills, with no 
matter fields, was analysed and a meromorphic /3-function was constructed which 
is finite at both weak and strong coupling. Up to a constant factor this /3-function 
reproduces the correct 1-loop Callan-Symanzik flow at both strong and weak coupling 
and interpolates between them analytically, although there is no a priori reason 
to interpret it as Callan-Symanzik /5-function away from the region of the fixed 
points. This analysis modified previous suggestions in the literature concerning the 
/3-functions for iV = 2 SUSY f|, § and evades the criticisms in ||. 

In the present paper the analysis is extended to include massless matter fields in 
the fundamental representation of SU(2) with Nf = 1, 2, 3 flavours. The construction 
uses the gauge invariant flow parameter u = tr < if 2 >, where tp is the Higgs field 
whose VEV is a free parameter, and the fact that the /3-functions are modular forms 
of weight —2 of a sub-group of the full modular group, T(l) w S7(2, Z)/Z 2 , depending 
on Nf. The significance of the parameter u was emphasised in where is was shown 
that u is the Legendre transform of the pre-potential. 
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Following m m a convenient choice of modular parameter for Nf > is 

9 8ni M _ 

T=~ + —, 1-1 

7T ^ 

where 9 is the usual topological parameter labelling #-vacua and g is the Yang-Mills 
coupling constant. In terms of r the relevant sub-groups of Sl(2, Z) for determining 
the /3-functions are: r°(2) for N f = 0; T(l) for N f = 1; r (2) for N f = 2; and r (4) 
for Nf = 3. 1 For Nf = and Nf = 2 these are larger than the monodromy group. 
The Nf = 1 case realises the full modular group, so the self-dual point r = i is a fixed 
point of the element sending r — > — i which is in the monodromy group. All the 
/5-functions discussed here do still have singularities somewhere in the fundamental 
domain of the relevant sub-group of r(l), they must do since any modular form 
of weight —2 must have at least one singularity within, or on the boundary of, 
the fundamental domain, but these singularities are off the real axis and, with the 
exception of Nf = 1, correspond to repulsive fixed points in both directions of the 
flow. For the Nf = 1 /3-functions there are two types of singularities off the real 
axis: one at r = i and its images under T(l) which is repulsive in both directions of 
the flow; and one at r = e i?r//3 and its images under T(l) which is attractive in the 
direction of decreasing Higgs VEV. 

The strategy is to use the technique of |2j where the Seiberg-Witten curves 
describing the various theories || are written both in terms of the 'bare' r = zoo 
and the renormalised finite r and the co-efficients compared to extract t(u). An 
important tool in the analysis is the modular symmetry derived in || |J where it was 
shown that N = 2 SUSY Yang-Mills has an infinite hierarchy of vacua with massless 
BPS states and the modular group relates these vacua to each other. The value of r 
in one vacuum is related to that of another by 

, , ar + b . 
where 7 = ( JGTc T(l) with T a sub-group of the full modular group 



c d / 

r(l) = PSl(2, Z). We therefore have, under a variation 5t of r, 



since ad — be — 1, so we expect that 

If (3 is meromorphic, it will be a modular form of weight -2 and this fact proves to 
be a powerful analytical tool. 

lr The notation is that of |Tc| ] and is summarised in the appendix for ease of reference. 
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In §2 the /3-function for Nf = discussed in |IJ is re-derived in terms of ( 1.1 



which differs from the normalisation in [|IJ. The cases Nf = 2, Nf = 3 and Nf = 1 
are treated in sections |3], ^ and |5| respectively where meromorphic /3-functions are 
proposed which vanish at all the strong coupling fixed points and are constant at the 
weak coupling fixed point. 

In principle the /3-functions for different values of Nf should be related by holo- 
morphic decoupling but the inclusion of non-zero masses for the matter multiplets 
makes analytic calculations much harder and cannot be pushed through using the 
techniques of the present analysis. In section 6 a perturbative approach is presented, 
using a strong coupling expansion and turning on a mass for one of the matter fields. 
It is shown that the limits td = — 1/t — > too and m — > do not commute, but never- 
theless a /3-function for the massive theory with acceptable behaviour near td = ioo 
can be constructed. 

Section 7 contains our conclusions. Two appendices give a summary of the 
conventions concerning Jacobi -^-functions and the technical aspects of the strong 
coupling instanton expansion used for the analysis in section 6. 



2. N f = 

This case was treated in QXJ using the normalisation appropriate to the adjoint repre- 
sentation of SU (2), f = ^ + ^r, but when matter in the fundamental representation 
of SU (2) is included it is better to define r = - + ||. In order to set the notation 
and illustrate the method for Nf ^ the derivation of the /3-function in |]J is given 
here using the original techniques of [fj|, adapted to the present notation. 
First recall the monodromy of the Nf = theory M. It is generated by 

M_\ !)• Mo 3 ™ 

with 

MocMo = M.2 1 where 

U (TWo leaves r = invariant, Ml 2 leaves r 
invariant). The two matrices .Moo and A4 , 

r^r + 4, r^——, (2.3) 
r + 1 

generate r°(4) and therefore j3 will be a modular form for r°(4) of weight -2. To 
determine its explicit form we start by following the method of 0. The massless 
Nf = 4 curve can be written 

y 2 = x 4 - 2uF(t)x 2 + u 2 , (2.4) 



M = J 0" (2 ' 2) 

= 2 invariant and .Moo leaves r = zoo 
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where 

F(T) _ m±^M (2 5) 

and the Nf = curve of || |9[ can be written as 

y 2 = x 4 - 2ux 2 + u 2 - A 4 . (2.6) 

Note that F(t) is an invariant function under T(2). 

Equation ( |2.4j ) is scale invariant in the sense that r does not depend on u. 
Equating co-efficients and eliminating u gives 

F ^ = -HT=T* (2 ' 7) 

hence 



u—F\t) = —F(F 2 - 1). (2.8) 



dr 
du 

where F' = Using (|A.ll ) to evaluate F'{j) from (2.5) gives 



P = 2m m3ii (2 9) 



u * = xm+m (2.io) 

du 7T 



finally leading to 



which is the result of 0, after taking account of the difference of notation (the 
variable r here is 2r in 0). The asymptotic behaviour of ( 2.10Q is 

u— — + — (2.11) 

which is the correct asymptotically free behaviour. However it was pointed out in 
|J that ( |2.10|) has a pathology in that it is singular at strong coupling, in particular 
at r = (u = Aq) and r = 2 (u = -Ag). 

A remedy was proposed in |l| , based on Q . The idea is to tame the singularities 
at strong coupling, without disturbing the asymptotic properties at u —>■ oo, by 
defining 

= Ju-Krw + KYJT 

v ' u m+n du v ' 

where m and n are positive integers (the sign is chosen so the direction of flow is the 
same as in 0). The correct behaviour at strong coupling then forces m = n = 1. 

The motivation behind ( 2.12|) relies on a theorem that any modular form of 
weight -2 for a sub-group T C T(l) can always be written as 

«r) = ^ (2.13) 
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where /(r) is a particular invariant function for T (essentially one with the smallest 
number of zeros plus number of poles, counting multiplicity) and P(f) and Q(f) are 
polynomials in / (see e.g [|Il|], page 111, the idea of applying this theorem to iV = 2 
SUSY was first proposed in |3j]). Now u(r) is just such an invariant function for r°(4) 
- the explicit form of u{r) follows from (|2.5 ) and ( |2.7| ) 



u _ ^(r)+^(r) 

Al 2^|(r)^(r) ' 1 ' } 



and it can be verified explicitly, using ( |A.5| ) and ( |A.6| ), that u is invariant under 
r°(4). Choosing the zeros of P and Q so as to get the correct asymptotic behaviour 
at r = 0, 2 and ioo without introducing unnecessary zeros or poles, one is led to ( [2.12j ) 
and, using m = n = 1 in conjunction with ( [2.10|) and (|2.14j ), this leads uniquely to 



|.-A|)(. + ^ = il 

u du TTi i/|(rj + i7|(r J 



The asymptotic behaviour is 





2 




-> — as 
iri 


r — 


/3(r)- 




as 




7^^ 


-2) 2 



(2.16) 



In particular the behaviour (3 ~ ^r 2 as r — > implies that /3(td) ~ - as To — >• «cxo, 
where = — 1/r is the dual coupling. Had any value of n been used in (|2.12j) other 
than unity, the result near r ?a would have been 

/3(td) oc (e 2 ^)"- 1 ! (2.17) 

which is not the correct asymptotic behaviour. An exactly parallel argument applies 
at r = 2 with n replaced by m. 

Near both the weak and the strong coupling fixed points this /3-function can be 
interpreted as a Callan-Symanzik /3-function for Nf = SUSY SU (2) Yang-Mills [|T| 
because u is a mass squared at weak coupling, where 

dr , 
f3 w — u— — , 2.18 
du 

while at strong coupling near r « 0, where w — Aq w AoOd with proportional to 
the mass of the BPS monopole, we have 

^-2 { u-M) d ^-2a D ^-. (2.19) 
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The flow is shown in figure 1. As the Higgs VEV is lowered there are attractive 
fixed points for all even integral r and repulsive fixed points for all odd integral r. 
There are repulsive fixed points in both flow directions at r = k + i for all odd k (/3(t) 
diverges at these points because i?|(t) = — i?§(r) there). From the transformation 
properties of the ^-functions in the appendix (|A.5|) and ( |A.6| ), the /3-function in 
(|2.15D is a modular form of weight -2 for the group generated by 



r^r + 2, t > r^T' ( 2 - 20 ) 
r + 1 

which is the group r°(2). r°(4) is a sub-group of r°(2) and the fact that the (3- 
function has a larger symmetry than that of the monodromy group is due to the 
Z2 symmetry of (3 under u — > —u, enforced by choosing m = n in (|2.12|) . This Z2 
is the anomaly-free residue of global i?-symmetry || which is not a symmetry of 
the effective action. Although u is not invariant under r°(2), it changes sign under 
t — > t + 2, u 2 is invariant. 

There are attractive fixed points at all the images of r = under r°(2) and 
repulsive fixed points at all the images of r = 1 under r°(2): at strong coupling all 
rational values r = q/m with q even are attractive fixed points and those with q odd 
are repulsive (with q and m mutually prime) as the Higgs VEV is decreased. 

The /^-functions in []IJ], using f = r/2 rather than r, are the same as ( |2.15| ) above 
as can be checked using ( |A.7|) from which 

i3(f) = 5« T > = h mrm - (22l) 

These are modular forms for T (2) which is generated by 

f->f + l, — - — , (2.22) 
2f + l' V ; 



equivalent to 



r^r + 2, r^^— (2.23) 

T + 1 



with r = 2f. The details of the flow generated by ( |2.21| ) were analysed in |T]|. 



3. Nf = 2 

The results of the last section can be immediately be used to guess the form of the 
/^-function for Nf = 2, which can then be verified explicitly. The monodromy for the 
Nf = 2 case is generated by 

M =( \ ?y m 00 =(~ i ~jy (3.1) 



2 1/' ™ V -1 
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with 

M oo M = M^ 1 where Mi=(~^ ^ (3.2) 

{M.q leaves r = invariant, M.\ leaves r = 1 invariant and M.^ leaves r = ioo 
invariant). The two matrices .Moo and -M^" 1 , 

T^T + 2, — - — , (3.3) 

2r + l' V ; 

generate T(2). If a /3-function is constructed which is invariant under the anomaly- 
free Z 2 acting on the w-plane, u — > —u, along the same lines as before then we expect 
it will have a further symmetry under 

r->r + l, (3.4) 

and so will be a modular form of To(2) of weight -2. Demanding the correct behaviour 
at t = ioo, r = and r = 1 leaves only one possibility and that is ( |2.21| ) with f 
replaced by r, namely 

2 1 

7rz i?|(r) + t?j(t) 

the normalisation being determined by the asymptotic condition 

(3{r) —> -. (3.6) 

We now verify ( |3.5| ) by explicit calculation using the same technique as in the 
previous section. The analysis initially parallels that of [f|: the massless Nj = 2 
curve is 

y 2 = x 4 - 2(u + 3A^/8)a; 2 + (u- A 2 2 /8) 2 . (3.7) 
Equating co-efficients with (|2.4|) gives 



or equivalent ly 



^ = (3J) 



Using ( |2.5| ) and ( |2.9| ), now leads to 



rfr (F-l)(F + 3) i flf(r)+fl«(r) 

du AF' 2vr 7?f(r)i?t(r) 1 ' 

as in 0. There are singularities at r = (where u = A|/8) and r = 1 (where 
u = — A|/8) which can be modified, along lines similar to §0, to give a /3-function 
with the correct strong coupling behaviour, 

{ u2 ~ dr 2 1 

P{T) ~ u 2 U Tu ~ «^(r)+^(r)' (3 - U) 

confirming ( |3.5| ). This flow is essentially the same as that of the Nf = 1 case treated 
in §0, except that r is rescaled by a factor of 1/2, and the flow shown in figure 2. 
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with 



4. N f = 3 

The monodromy for Nf = 3 is generated by 

A<o-Moo = -Mli/2 where ^-1/2 = ^1 (4-2) 

(A^o leaves r = invariant, leaves r = — 1/2 invariant and .Moo leaves r = zoo 

invariant), [[J. The two matrices .Moo and M-q 1 , 

r-r + 1, r^^, (4.3) 

generate To (4). But it is not obvious what the full symmetry of the /3-functions 
might be as Z2 does not play any role in the Nf = 3 theory ||. We must therefore 
perform the explicit calculation. Our starting point this time is the massless Nf = 4 
curve in the original form of ||, namely 

y 2 = x 3 - ^g 2 {r)u 2 x - ^g 3 {r)u 3 (4.4) 



with 



and 



92(r) = l(^(r) + ^(T) + ^(r)) (4.5) 



9*(t) = - (*J(r) + <(r)) (^(r) + 0*(r)) (<(r) - 0*( T )) , (4.6) 
the combination 

= JCr) (4.7) 



^ 3 (r) - 27^(r) 54^^f 

being Klein's absolute invariant which is invariant under the action of T(l) on r. 
The massless Nf = 3 curve is 

V 2 = (x 2 -^(x-u)^j (x-u). (4.8) 

Without loss of generality we can set 2 A 2 . = 64 and then eliminate the quadratic 
term in (|4.8j ) by shifting x ^ x + (u + l)/3 giving 

y 2 = x 3 — \u 2 -Au + l)x - 4=(2m - l){u 2 + 8u - 2). (4.9) 
3 27 



2 This is equivalent to defining the dimensionless variable u — 64u/(A§)) and then dropping the 
tilde. 
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There are singular points were the roots of this equation coincide, at w = 0, it = 1/4 
and \u\ = oo, corresponding to r = 0, r = — 1/2 and r = ioo respectively. 
Equating co-efficients between ( [4.4[ ) and ( |4.9| ) and eliminating u yields 

4 (u 2 - Au + l) 3 , 
J « = 27 \,*(1-4 M ) - < 4 ' 10) 

To extract the /3-function from this equation we will need an explicit expression 
for u{t). To this end define 

in terms of which (f4.10| ) reads 

(i + y) 3 _ (i -x + x 2 f 

V 2 ~ (1 -X) 2 X 2 

with the three roots 



(4.12) 



y = -X(l-X), Y 2 = { ±^P- and y 3 = , (4.13) 
Comparing with the three asymptotic forms 

• oo, Y — >• oo, 

3, y - 0, (4.14) 

r->-l/2, X->1, u->l/4, y^oo, 

only Y3 has the correct asymptotic behaviour at the three singular points and so we 
must choose 

u 2 X -&%d\ 



— > ioo, 


X - 


+ 1, 


\u 




X - 




u 


-1/2, 


X - 


-1, 


u 



1-Au (1-X) 2 " 0f ' (4 ' 15) 
which is an invariant function for T (2) [[!]]]. Solving for u the asymptotic conditions 
pick out the unique solution 

which is, of course, an invariant function for T (4). 

Differentiating this equation with respect to r, and employing ( |A.11|) , yields 

dr 2i 1 , 
u— = x. 4.17 

du 7r(^|(r)+^(r)) 2 

This has the correct asymptotic form asr-> ioo, 



and is well behaved at r = 

dr 2i , / , s 

u— — > r 2 , 4.19 

but diverges at r = —1/2. 

The method used in the previous sections for Nf = and Nf = 2 is not im- 
mediately applicable here since one of the fixed points is at u = 0: to eliminate the 
singularity at u = 1/4 we would have to multiply by (u— l/4) m /u m for some positive 
m which would introduce another singularity at u — 0. One strategy is to use m — 1 
to remove the singularity at u = 1/4 and shift the other fixed point away from the 
origin by shifting u: thus let v! = u — e for some constant e and define 

= >' + f -l/4)K + e)rfl (4 . 20) 

This does not disturb the behaviour at \u\ — oo. Clearly this is equivalent to 

. , (u- 1/4) dr 

/? r = - y - Liu— 4.21 

[u — e) du 

which preserves the good behaviour at r = ioo and r = and gives the correct 
behaviour at r = —1/2. 
In terms of -d-functions 



P(T) = — ; — r — ; , — t ; n0/ xx2 , (4-22) 



1 

2ttz i^(r)^(r) + e^|(r) - ^( r )) 
and has the following asymptotic forms 

1 



/3{t) — > (4.23) 
2 (r + 1/2) 2 



-i/2 iri (1 — 4e) 



If e is real and < e < 1/4, the /5-function is finite at both r = and r = —1/2 
and flows in the right direction, i.e. in towards the fixed points as the Higgs VEV is 
lowered. 

In fact just such a constant shift of u was found to be necessary in instanton 



calculations performed to check the validity of the Seiberg-Witten curve |JL2| , [13 
These instanton calculations give e = 4/27 when A 3 /8 is set to one as we have done 
here. The point is that, since Nf = 3 has no discrete symmetry acting on the u- 
plane, there is no a priori way, using the techniques in 0, to determine where the 
origin of the u-plane should lie and a constant finite shift of u does not affect the 
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weak coupling physics. Seiberg and Witten chose u = to correspond to r = but 
the instanton calculations show that u' = —4/27 is a more natural choice of origin. 

The flow (|4.22|) with e = 4/27 is plotted in figure 3. At first glance it looks very 
like the flow for Nj = and Nf = 2, with r rescaled, but closer examination reveals 
subtle differences. Figure 3 is not symmetric under r — > r+1/2, it is slightly distorted 
and the repulsive fixed point close to r = (1 + i)/4 is not exactly at r = (1 + z)/4, 
it is displaced away from the top of the semi-circular arch by a small amount. The 
flow looks like a distorted version of the Nf = 2 flow with r rescaled by a factor of 2, 
T N f =2 = 2tjv /=3 . This is because r°(4) acting on r is equivalent to T(2) acting on 2r: 
the Nf = 3 /3-functions for 2r are therefore modular forms of T(2). Modular forms 
of T(2) can be obtained by distorting modular forms of r (2). Provided < e < 1/4 
the unstable fixed point of the Nf = 3 flow lies on the semi-circle in the upper-half 
r-plane spanning the two points r = and r — 1/2 on the real axis. The special case 
e = 1/8 has a higher symmetry than other values because this corresponds to the 
unstable fixed point being at the top of the semi-circular arch at r = (1 + i)/4 and 
this gives To (2) symmetry acting on 2r. Other values of e have the lower symmetry 
of T(2), as in figure 3. Flow diagrams like figure 3 have been postulated for the 
quantum Hall effect when the electron spins are poorly split |L4 . 



5. N f = 1 

The Nf = 1 case has been left till last because it is more involved than the three 
cases already considered, though paradoxically it has a higher symmetry - - the 
monodromy generates the full modular group r(l). The monodromy is calculated 
in |J: there are four singular points, at r = 0, r = 1, r = 2 and r = zoo, with 
monodromies 

M =(\ ?V Ml '° 1 



117' V -1 2 



respectively. Now 



and 



M >=[-i 3)' ^("o 1 IiJ (") 



MlM 1 =l® 1 Jj (5.2) 



< 1 1 



MoMxMo 1 = ( ! J (5-3) 



so the two operations 



t->-1/t and r^r+1 (5.4) 
are in the group generated by ( |5.1| ), which is therefore the full modular group T(l). 
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Again the discussion starts along the lines of @. We take the massless Nf = 4 
curve 

y 2 = x 3 - -g 2 (r)u 2 x - -# 3 (^ 3 (5.5) 

with 

g i(r) _ (n+n+nr _ 7(r) 

gi(r)-27gi(r) 54^8 W 
and compare this with the massless Nf = 1 curve 

y 2 = x 2 (x-u) - A^/64. (5.7) 

First eliminate the x 2 term in the Nf = 1 curve by shifting x — > x + u/3, 

y 2 = x 3 x-i —u 3 + —). 5.8 

y 3 V27 64/ v ; 

Equating co-efficients with ( |5.5| ) and eliminating u then produces 

J(r) = -\-^TT, (5-9) 
4 zr + 1 

where the scale has been set by choosing 27A^/256 = 1. As r — > zoo, J — > oo and 
this corresponds to w — > — oo. But J is invariant under T(l) and so J —>■ oo at the 
three points r = 0, r = 1 and r = 2 as well and these correspond to the three roots 
of u 3 = — 1 . Differentiating equation (5.9) then leads to 

^ = 3 l#nJj 7 (5 - 10) 



and QA.llj) gives 



As r — > zoo, 



J f ■&% + 0f + &l 



J' 2m (^ + tft) (0* - z?t) ($j + z?|) ' 



(5.11) 



J/J'-f and u^ = ^, (5.12) 
27r du 2tc 

which is the correct asymptotic behaviour. But J/ J' has the same value at r = 0, 1 

and 2 as at r = zoo and u 3 = — 1 at these 3 points, so (|5.10|) diverges at strong 

coupling. Following the same procedure as before the singularities in ( 5.10Q can be 

eliminated, without disturbing the behaviour at u ~ — oo, by using 

u 3 + 1\ dr f u 3 + 2\ J 



This is not a modular form for T(l) however, since u is not invariant. Solving (JIT 
for u gives 



u 3 = -2 (j±y/j(J-l) ) . (5.14) 
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u 



and both roots are necessary: the upper sign for u — > — oo and the lower sign for 
-1. Eliminating u 3 from (|5.13|) then gives an ambiguity in the direction of flow 



u 3 J du J' i 7l ^2(dl + dl + §\) 

The resolution of this problem is that the true /3-function, which should be 
a modular form for must have poles or zeros that are not accounted for in 

( |5.13|) . This equation was derived by making the minimal modification of ( 5.10 ) that 
would remove the infinities at m 3 = —1 yet not disturb the asymptotic behaviour 
at u — — oo. For Nf = 1,2 and 3 this minimalist approach worked. For Nf = 1 
it does not, there must be another pole or zero somewhere. In fact ( |5.15| ) already 
has a singularity at r = e t7T ^ 3 , where $| + $3 + $4 = corresponding to u = 0, and 
its images. This is a fixed point of T(l), meaning that there is at least one element 
7 G T(l) that leaves it invariant. 3 The point r — i, where d\ = $f, is also a fixed 
point of T(l) and, if the flow commutes with T(l), then t — i should be a fixed point 
of the flow also. The next simplest assumption that can be made is that the only 
fixed points above the real axis are at r = i and r = e m l 3 and their images under 
r(l), and we shall therefore assume that these are the only places above the real axis 
where a pole or a zero of the /3-function is allowed. We shall further assume that 
the total number of poles plus zeros (counting multiplicity) is the smallest number 
compatible with these requirements, since this was our experience for Nf = 0, 1 and 
3 and was an assumption in ||. Jacobi's invariant function takes the values J = at 
r = e m l 3 and J = 1 at r = i. Therefore, with the above assumption, the /3-function 
must be of the form 

P(r) = -3 J {J r 1} (5.16) 

with m and n integers. In order that the /^-function has the correct asymptotic form 
as r — > ioo it must be the case that m + n = 1 and minimising the total number 
of zeros and poles leaves only two possibilities m = 0, n = 1 or m = 1, n = 0. 
Examining these two possibilities the first gives 

j_l _ 2 ^l + ^l-^l + ^j) 
J' iri (i + i)| + * 2 1 ' 



'2 ~ ^3 ~ U 4J 

which is singular when + $3 + 1?| = 0, that is at r = e ln ^ 3 and its images, and the 
second gives 

J' 27Ti(^| + ^)(^-^)(^ + ^) ^ 

which is singular when (i?| + i?!)^! — ^4) (^3 + ^4) = O5 th.a.t is at r = i and its 
images. 4 The latter has a milder singularity, and clearly has a smaller number of 

3 Actually there two such elements: 7 = f ^ ^ J and 7 = ^ ^ ^ 
H\ = §\ at r = i, <&% = -d\ at t = (1 + i)/2 and = -<&% at r = 1 + %■ 
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poles plus zeros than the former, and so is the unique choice that fits the criteria. 
So we conjecture that the analytic /^-function for Nf = 1 is 



3 



0(t) = - ( ^\ = 3 (g|±j§±jg) (5 19) 

which differs from ( 5.13Q only by using u 3 +2 in the denominator, rather than u 3 . This 
flow is plotted in figure 4, the pole at r = i, corresponding to u 3 = —2, renders this a 
repulsive fixed point in both flow directions while the zero at r = e tn ^ 3 , corresponding 
to u = 0, is an attractive fixed point in the direction of decreasing Higgs VEV. 



6. Massive Matter Multiplets 

When the matter fields in the fundamental representation of SU(2) have non-zero 
mass the analysis is much harder. For massive matter multiplets the Seiberg-Witten 
curves determine r(u, rrii, Ajvy) with i = 1, . . . , Nf and the Nf — 1 case can be deter- 
mined from the Nf case by holomorphic decoupling ||: set = m Nf A 4 Nf Nf and 
send rriNf — > oo and A^ f — > keeping A^-i finite. In principle therefore it ought 
to be possible to distort figure 3 for example, by turning on one mass, and recover 
figure 2 as the mass goes to infinity. This would correspond to a family of /^-functions 
obtained by differentiating t(u, m^A N ) with respect to u and considering the to 
parameterise different /5-functions: near the fixed points this corresponds to defining 
/^-functions by varying the W /± -boson, and therefore also the gluino, masses while 
keeping the quark masses fixed. 

In this section we address the question of quark masses using a strong coupling 
expansion, taking the Nf = 3 case for illustrative purposes. The details of the 
analysis are rather technical and so are relegated to appendix B, from which we 
quote the relevant formulae here. 

For Nf = 3 with three different masses mi, m2 and 777.3 the Seiberg-Witten curve 
is uniquely determined by 777.1, 777 2 , 7773 and the r (4) invariant parameter u. With 
the term quadratic in x eliminated the curve is 

y 2 = X 3 — — (iL 2 — AuAl + A3 + 3(777^ + 7772 + m l)^-3 ~ 6A3777! 77727773 j X 

- 1 {(2u - Al)(u 2 + 8u~A 2 3 - 2Al - 9(777? + 777^ + m^) (6.1) 

— 18(w + A3)A3777i77727773 + 27A\(m\m\ + m\m\ + 7773777?) | , 

where, in the notation of U, A 3 = A3/8. In principle this curve determines r{u,mi) 
and there are four different masses that could be varied to define /5-functions. To 
connect with the /3-functions of §|4] we shall keep rrii fixed and vary only u. 

In the massless case the strategy was to find an explicit expression for u(r), 
in terms of Jacobi 79-functions, and then use the properties of the "^-functions to 
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determine The problem in the massive case is to find u(t, rrii) and this is much 
harder. We shall first simplify the problem and only consider one non-zero mass 
mi — m i m 2 — m 3 — 0. Even then one cannot hope for a closed form solution. 
However it is shown in appendix |D| that, at strong coupling where T£> = — 1/t — > 200, 
it is appropriate to expand in qn := e i7TTD ^ 2 . The details are somewhat technical and 
left to the appendix but for one non-zero mass there are singularities at r = when 
u = ±mA 3 . Near u = mA 3 the expansion 



u 



(td) = mA 3 (1 + a(m)e iwTD + ■ ■ ■) (6.2) 



is derived in appendix B where a(m) is an unknown function of m/A 3 which cannot 
be determined without further assumptions, but it diverges like ~ 1/m 2 as m — > 
and a(m) — > 16 as m — > oo. 

For u infinitesimally close to m the /3-function will be of the form 

Or 

P(t) oc -(u-m)— , (6.3) 
ou 

assuming the BPS monopoles have mass oc u — m. Using 

du 



dr D 

gives, for u near m, 



%TimK^a{m)e mTD ~ m(u — m), (6-4) 



P(t d ) oc -(u-m)^- « -, (6.5) 

OU 7T 

which is perfectly well behaved, even in the massless limit m —> 0. Indeed 

IT 2 

(3{t) oc — (6.6) 



7T 



which, up to a constant, is the same behaviour as equation (|4.23|) even though it is 
clear from ( |6.2j ) that the limits td —> ioo and m — > do not commute, since a(m) 
behaves as 1/m 2 as m — > 0. 



7. Conclusions 

Explicit expressions have been proposed for the /3-functions of = 2 SUSY SU (2) 
Yang-Mills with massless matter fields in the fundamental representation. Asymp- 
totically close to the strong and weak coupling fixed points they co-incide with the 
1-loop Callan-Symanzik /3-functions, up to a constant factor. 

The /^-functions are modular forms of sub-groups of T(l) for each value of Nf. 

N f = r°(2) 

N f = i r(i) 
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N f = 2 r (2) 
N f = 3 r (4), 



(7.1) 



for Nf = 1 and Nf = 3 the group is the same as the monodromy group and for 
Nf = and Nf = 2 it is larger, due to the Z2 action on the w-plane. 

The /5-functions are determined by demanding that they have the correct asymp- 
totic behaviour at both weak and strong coupling fixed points. The relevant flows, in 
the direction of decreasing Higgs VEV, are shown in figures 1,4,2 and 3 respectively. 

These functions differ from previous expressions in the literature in that they 
have the correct behaviour at all the strong coupling fixed points. In all cases the 
/5-functions have a singularity in the interior, or on the boundary, of the fundamental 
domain, corresponding to an unstable fixed point which is repulsive in both directions 
of flow. This is a necessary consequence of their being modular forms of weight -2. 
For the Nf = 1 case, for which the monodromy group is the full modular group this 
repulsive fixed point is at the self-dual point r — i and its images and in this case 
there is also a fully attractive fixed point at r = e ln ^ 3 and its images. The physical 
significance, if any, of this attractive fixed point is not clear. 

The case of finite masses is probably intractable using the methods developed 
here, though it may well be possible to make progress with other gauge groups by 
using the methods in ||. Nevertheless it has been possible to show, in a strong 
coupling expansion, that turning on one quark mass in the Nf = 3 case still allows 
for a /3-function with the correct asymptotic behaviour near r = 0. Unfortunately it 
is not possible to say anything about the /3-function away from r = because of the 
limitations of the technique. 

It is a pleasure to thank the Perimeter Institute, Waterloo, where this work was 
completed, for hospitality. This work was supported in part by Enterprise Ireland 
Basic Research Grant no. SC/2003/415. 



A. Appendix: properties of Jacobi ^-functions 

We collect together some useful properties of ^-functions. The definitions are those 
of [|IIJ and most of the formulae here are proven in that reference. The three Jacobi 
^-functions used in the text are defined as 



00 



# 2 ( r ) = 2^g^ 2 =2gi J](l-g 2n )(l + g 2 ") 2 , (A.l) 

n=0 n=l 

00 00 

Mr) = J2 = II^ 1 - ^X 1 + ^) 2 ' ( A - 2 ) 

n=— 00 n=l 

00 00 

Mr) = E ("W = II( 1 - ^ 2 ")( 1 - ( A - 3 ) 



n=l 
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where q := e l7TT . 

These three ^-functions are not independent but are related by 

The following relations can be used to determine their properties under modular 
t r ansf ormat ions : 

Mr + 1) = e^ 2 (r), tf 3 (r + 1) = Mr), Mr + 1) = Mr), (A.5) 



# 2 (-1/t) = V^Mr), 

M-i/r) = V^Mr), (A.6) 
M-^/r) = V^Mr)- 

The duplication formulae 

^(2r) = i((^(r)-^(r)), 

$i(2r) = \(Mr) + Mr)) , (A.7) 
^(2r) = ^(t)^(t), 

and 

tf 3 (4r) = -(tfsW + ^r)), tf 2 (4r) = - (<? 3 (t) - 4 (t)) (A.8) 
are also useful. 

At the special points r = e* ir ' 2 and r = e" 1 "/ 3 the ^-functions have the values 

Me ln/2 ) = V2Me l7T/2 ) = V2Me ln/2 ) = -K (sin (|)) , (A.9) 

e- ff / 4 ^(e" /3 ) = e" l7r/12 ^( e -/ 3 ) = ^-&\{e^) = (sin (^)) , 

where i^(A;) is the complete elliptic of the second kind: K (sin(ir / 4)) = ^^(r(l/4)) 2 , 
with T(l/4) rs 3.6256 the Euler T-function evaluated at 1/4, and X(sin(7r/12)) ~ 
1.5981. 

The ^-functions have the following asymptotic forms 

r -> ioo : i? 2 (t) « 2 e^ 1 -> 0, 3 (r) -> 1, $ 4 (t) -> 1; (A.10) 

r^O : tf 2 (r)«J? tf 3 (r) « JI, tf 4 (r) « 2^ e"£ - 0. 



In addition they satisfy the following differential equations (see [fLl] , p. 231, equation 
(7.2.17)), 

1?', ^ Z7T .„ 



d b 

In the text T (N) consists of matrices 7 = ( ) in m PSl(2, Z) with 



c = mod N, sometimes written 



c d 



* 

T°(N) consists of matrices with 6 = mod N 



7=1* * ) modiV, (A.12) 



7=1* ° ) mod N (A.13) 



and T(N) consists of matrices with 

7= ( J mod N. (A.14) 

B. Appendix: Massive Matter Fields for Nf = 3 

The starting point is the curve for four massive matter multiplets in the fundamental 
representation ||, 

y 2 = (x 2 — c\u 2 )(x — Ciu) — c 2 (x — Ciu) 2 A — c 2 (c 2 — c%)(x — C\u)B 

+ 2c 2 (c 2 - c 2 2 )( Cl x - c 2 «)C< - c 2 (c 2 - c 2 ) 2 A (B.l) 

with 

ci(r) = -(0*(r) + tft(r)), c 2 (r) = -(0«(r) - <5>*(r)), (B.2) 

and 

1 



A = ^^m 2 , B = ^^m 2 m 2 , C = m 1 m2m 3 m4 D= r//~ iir-mj.. 



- 2 - 2 - 2 

i i<jr i<j<k 

(B.3) 

First eliminate the quadratic term in x by shifting x — ■> x + (ciW + C2A)/3, giving 
y 2 = x 3 - -P(r, u, m,)x - u, mj), (B.4) 

where 

P{T,u,fhi) = (c 2 + 3c 2 )m 2 -4c 1 c 2 2 Au + 3c 2 {c 2 1 -c 2 2 ){c 2 B-2c 1 C) + 4A 2 , (B.5) 

is quadratic in u and 

Q(r,u,mi) = 2 Cl {cl-9c 2 2 )u 3 + 3c 2 2 {5cl + 3c 2 2 )Au 2 (B.6) 
- 2 (6cic^ 2 + 9cic^(c 2 - c\)B + 9c 2 (c 2 - c^)(c 2 - 3c 2 2 )C) u 
+ 2c^ 3 + 27c 2 (c 2 - c\) 2 D + 9c 3 (c 2 - c 2 2 ){c 2 B - 2c 1 C)A 
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is cubic in u. 

Next this curve can be reduced to the Nf = 3 curve using the holomorphic 
decoupling of ||: send m 4 — > oo and take the 'bare' coupling r — ► ioo, so c 2 — ► 
and C\ — > 1, keeping C2m 4 = A 3 /8 := A 3 finite. Taking this limit, and dropping the 
bars on u, m\, m 2 and 777.3, gives the massive Afy = 3 curve, 

y 2 = x 3 - -(u 2 - 4wA 2 + A3 + 3BA 2 3 - 6CA 3 )x (B.7) 
3 

- ^(2u - A 2 3 )(u 2 + 8 M A 2 - 2A 4 3 - 95A 2 ) + (it + A^)CA 3 - D, 

with 

I? = m 2 + m\ + m 3 , (7 = m 1 m 2 m 3 , D = m\m\ + m 2 m 2 + m^m\. (B.8) 

This reduces to equation ( |4.9[) when the three masses are set to zero and A 3 = 1. 
Following this curve is now compared to (B^) with m 4 = 0, u independent of r, 
that is equation ( [B.4j) - (|B.6|) with (7 = and 

A = rh\ + fn\ + m 3 , S = m\m\ -Hm 2 ,???,^ +m^fhf, D = m\m\m\. (B.9) 

Equating co-efficients leads to 

(c^ + Sc^^-^ic^m + Sc^c^-c^E + c^ 2 = it 2 -4iiA2+A^+3M 2 -6CA 3 (B.10) 



2ci(c 2 - 9c^)u 3 + 3c|(5c 2 + 3c 2 2 )Au 2 - 2 (6dc 4 2 A 2 + §c x c\{c\ - c\)B) u 

+2c 6 2 A 3 + 27c 2 2 (c\ - c\) 2 D + 94(cj - c\)BA (B.ll) 
= (2u - A 2 )(w 2 + 8wA 2 - 2A\ - 95A 2 ) - 18(u + A 2 )<7A 3 + 27 D. 

In the massless case we can eliminate u from these two equations and determine u(t), 
but now we want to determine u(r, m,) with w, A, -B and D all unknown functions 
so there is not enough information to solve the problem completely. Nevertheless 
we can still get information from these equations. The symmetry group for Nf = 3 
is r (4) and u, rrii and A 3 should be invariants of T (4) so the right hand sides of 
( |B.10| ) and ( |B.11| ) are also invariants. Using the transformation properties of the 
^-functions, ( |A.5| ) and ( |A.7| ), this implies that u, A, B, and D are modular forms of 
weights -2, -4, -8 and -12 respectively, so fhi have weight -2. 

For simplicity we shall focus on the case of a single mass, m\ = m, m 2 = m 3 = 

so 

A = m 2 , B = 0, D = (B.12) 
B = m 2 , D = 0, (7 = (B.13) 
and the equations simplify to 

(c 2 + 3c 2 )u 2 - A Cl c 2 2 m 2 u + c\m A = u 2 - 4uA 2 + A3 + 3m 2 A 2 (B.14) 
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2ci(4 - 9c l)u 3 + 3c 2 2 (5cl + 3c 2 2 )m 2 u 2 - 12 Cl 4m 4 u + 1c\m % 

= (2m - A 2 )(w 2 + 8uA 2 3 - 2A 4 3 - 9m 2 A3). (B.15) 

We now want to eliminate u and rh to get u(r,m) but there is still not enough 
information. However, knowing that u and rh 2 are modular forms of weight -2 and 
-4 respectively, we can say something about their functional form at strong coupling. 
To see how this works let us first look at the case rh = 0, where the explicit solution 
is given in §f|. Using the details there one finds 

A 2 

u(t) = 5 (B.16) 

(^l(T)-^(r)) 2 

which is indeed a modular form for To (4) of weight -2 and it vanishes at strong 
coupling, t — > 0. Now r ^ — 1/t = td is not in To (4), so u is not a modular form 
under this transformation, rather 

T D (^i(TD) -^I(td)) 

Writing = e* 711 "- / 2 we have the strong coupling expansions 

A 2 

u{t d ) = -§■ (1 + 8q D + 40? D + 160f D + ■•■). (B.18) 



T 



D 



and, from Q4.16Q , 

<r D ) = - ( J; {T ^ D \, 2 = (1 + 8q D + U? D + + •■■)• (B.19) 

For non-zero m it is consistent with all we know to assume a similar form 

A 2 



(td) = Jr{u + u x q D + u 2 q 2 D + u 3 f D + •••). (B.20) 

and similarly for rh 2 



A 2 

•• 2 'r D ) = (o + ai<?D + a2«?D + 03<7d H ) > (B.21) 



rn 



T 



D 



where Uk and at are functions of 777./A3, with vanishing for m = 0. A similar 
strong coupling expansion for u has no prefactor of 1/t|, because w has weight zero 
not -2, 

— = UQ J r Ul q D + U2 g2 + U3 ^3 ^ _ (B.22) 

Using these expansions in ( |B.14| ) and ( |B.15| ), together with 

Cl (r) = Cl (-l/r D ) = -$($ 3 (t d ) 4 + Utd) 4 ), 



C2 (r) = C2 (-l/r D ) = - T 1o,(t d )\ (B.23) 
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we can equate powers of qo to obtain recurrence relations between the Uk and the dk- 
Without making further assumptions there is not enough information to determine 
u(t, m) but we can still extract useful information about the strong coupling (3- 
function. At zeroth order in q^ ( B.14|) and ( |B.15| ) (with A 3 = 1 for simplicity) yield 
three possibilities: 

u = ±m, uo = 1 =F 2m — — ; and Uq = m + -,u = m — . (B.24) 

Only the first two are relevant for the strong coupling fixed point, r = at u = 0, of 
the massless Nf = 3 theory (in the massless case u — 1/4 is associated with r = —1/2 
where td — 2). Choosing the root u — m, at order q^ we find a pair of equations 
linear in u\ and u\ which are degenerate. Solving for u\ gives 

fix ;= _ (;»- 2 )"i (B.25) 
1 (2m -1) 4 v 7 

At order the pair of linear equations for for u 2 and U2 are parallel in the u<i — u 2 
plane and have no solution unless they co-incide, which only happens if 

= 0. (B.26) 
2m - 1 y J 

For m = this is automatic, but for m ^ it forces U\ = 0. Setting u\ = then 
gives one linear equation relating m 2 to m 2 , 

2(m - 2)u 2 + a (3a + 8m - 4) a 2 /r> n _v 

U2 = 2(2^ T (R27) 

At order g^, the two linear equations for u% and u% are again degenerate giving only 
one constraint which can be used to solve for U3 

(m-2)w 3 + ai(3a + 4m-2) a 3 /D no , 

U3 = (2^i) T (B - 28) 

At order q\y one again obtains two parallel lines in the 1*4 — plane which do not 
intersect unless 

(d + 8m - 4) 2 , , 

u 2 = B.29 

2m 

in which case U4 can be obtained as a function of U4, d , a 2 and m. 

One can continue but for the present purposes we have gone as far as necessary. 
We are only really interested in u and we have 

, . (do + 8m — 4) 2 n 

u{t d ) = m ± } -q 2 D + ■■■ B.30 

2m 

with do and undetermined function of m, but independent of td- To get the dimen- 
sions correct we should re-instate A3 and write 



u 



y D ) = m A 3 (1 + a (m)e i7ITD + • • •) (B.31) 
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where 



a{m) := ± (M§ + 8mA, -4Ag)' 



Note that, in order to pin down the co-efficient ui one has to go to order q% and 
examine U4. At every value of k in the expansion one gets a pair of linear equations 
in Uk and in terms of m and the a^' with k' < k. For k = 1 these equations are 
degenerate and Mi is not determined; for k = 2 the equations have no solution unless 
iii = in which case they are again degenerate and u 2 is undetermined; for k = 3 the 
equations are again degenerate and % is undetermined and for k = 4 the equations 
have no solution unless u 2 has one of the two possible values shown in (|B.31[) . 



The explicit from of a(m) is not needed in the analysis, but we can fix its 
asymptotic form as m — > 00 using holomorphic decoupling ||. For Nf = 2 equation 
JD gives 



and this should agree with ( |B.31[ ) as m — > 00, A 3 — > with A 2 , = mA 3 = 8mA 3 fixed. 



Hence a(m) — > 16 as m — > 00. 

A 2 

The other singularity of the Nf = 2 theory, at u = is obtained by holo- 

morphic decoupling in the weak coupling limit of the Nf = 3 theory by expanding 
around u = —m. 

Notice that the m = expansion for u in equation (|B.19|) contains a term linear 
in q D while ( |B.31| ) does not. This is because, when m 7^ 0, ( |B.26[ ) forces us to set 



Ui = 0. Since do — > as m — > 0, a(m) diverges like SAg/m 2 as m — > and the limits 
7~d — > ioo and m — > do not commute. 

We can perform a similar expansion around r = —1/2, using w = m 2 + 1/4, 
where td = 2 + is with e small. The analysis is simpler than the uq = ±m case in 
that at each order, at least up to order 5 in qjj which is as far as we have gone, one 
simply finds a pair of linear equations in and which can be solved in terms of 
dk' with k' < k. The details are omitted but one finds 

(a -4m 2 + l) 4 a 1 (a -4m 2 + l) 3 

Ul = u 2 = u 3 = 0, u A = -4 — — — — — , u 5 = -lQ — — — , 

{4m z — l) z (4m z — l) z 

(B.34) 

and so, for m = 2 + is with e small, 

M(TD ) = m' + i-4 (a " Y ^ , -- 2 ''+^ M T^ m \t 1)3 e-^+- ■ ■ . (B.35) 
v 1 4 (4m 2 - l) 2 (4m 2 - l) 2 ' 

Re-instating A 3 now gives 

u{t d ) = m 2 + ^ + A 2 a(m)e- 27r£ • • • . (B.36) 
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where 

, . (a As - 4m 2 + A 2 ) 4 

a(m) = -4^^ - 3J . B.37 

V ; (4m 2 -A 2 ) 2 Af v ; 

As m — > oo this point in the w-plane goes out to infinity in the Nf = 2 theory. 

Equations ( |B.31|) and ( |B.36|) were the aim of this appendix and are used in in 
the discussion of the /5-functions for the massive Nf = 3 theory at strong coupling. 
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Fig. 1: Flow of effective coupling of iV = 2 SUSY Yang-Mills with N f — 0. The 
arrows indicate the direction of the flow as as the Higgs VEV is reduced. The (3- 
functions are modular forms of r°(2) and the pattern repeats under r — > r + 2. 
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Fig. 2: Flow of effective coupling of iV = 2 SUSY Yang-Mills with N f = 2. The 
arrows indicate the direction of the flow as as the Higgs VEV is reduced. The (3- 
functions are modular forms of r (2) and the pattern repeats under r — > r + 1. 
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Fig. 3: Flow of effective coupling of iV = 2 SUSY Yang-Mills with N f = 3 and 
e = 4/27. The arrows indicate the direction of the flow as as the Higgs VEV is 
reduced. The /^-functions are modular forms of r (4) and the pattern repeats under 

T -»• T + 1. 
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Fig. 4: Flow of effective coupling of iV = 2 SUSY Yang-Mills with Nf — 1. The 
arrows indicate the direction of the flow as as the Higgs VEV is reduced. The (3- 
functions are modular forms of r(l) and the pattern repeats under r — > r + 1. 
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